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1 Introduction 

Inflation produces a vast ensemble of infrared scalars and gravitons. This is thought to be the source 
of primordial scalar and tensor perturbations [1]. It is natural to wonder what effect these ensembles 
have on other particles. That sort of question can be answered by computing the scalar or graviton 
contribution to the appropriate IPI (one-particle-irreducible) 2-point function and then using that 
result to quantum-correct the linearized field equation for the particle in question. 

The IPI 2-point function for a scalar is known as its "self-mass-squared", —iM'^{x;x'), and the 
quantum-corrected, linearized field equation for a massless, minimally coupled scalar is. 
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where g^u{x) is the spacelike metric tensor. The fermion's IPI 2-point function is called its "self- 
energy", ~i[iTij]{x;x'), and the quantum-corrected, linearized field equation for a massless fermion 
is, 
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where e^a(x) is the vierbein field, 7j^- are the gamma matrices, A^hc{x) is the spin connection, and 
J^'^ = — 1[7^,7'^] are the Lorentz generators. The IPI 2-point function for a photon has the evocative 
name "vacuum polarization", +i[^Il'^]{x; x') , and the quantum-corrected Maxwell equation is. 
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where F^i, = d^A^ — di^A^ is the field strength tensor and J^{x) is the current density. And the 
IPI 2-point function for a graviton is termed the "graviton self-energy", —i[^'^Ti^'^]{x;x'), and the 
quantum-corrected, linearized Einstein equation is. 
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where V'^"'^" is the Lichnerowicz operator, k^ = IGttG is the loop counting parameter of quantum 
gravity and T^^^ is the linearized stress tensor. 

Many results of this type have been derived in recent years, with the background geometry of 
primordial inflation modeled using the cosmological patch of de Sitter space. The effects of massless, 
minimally coupled (MMC) scalars are simplest to study. A quartic self-interaction leads MMC scalars 
to develop a growing mass [2]. The vacuum polarization from charged MMC scalars causes the photon 
to develop a mass [3] and engenders profound changes in electrodynamic forces [4]. MMC scalars which 
are Yukawa-coupled to a fermion make the fermion develop a growing mass [5]. And MMC scalars 
do not have any effect on gravitons provided one can absorb certain surface terms into perturbative 
corrections of the initial state [6] . 

The effects of inflationary gravitons are more difficult to work out (everything is tougher in quan- 
tum gravity!), but studies have been made of what they do to fermions and to MMC scalars. The 
results are interestingly different: whereas inflationary gravitons induce a slow growth of the fermion 
field strength [7] they have no secular effect on MMC scalars [8]. The difference seems to be due 
to spin. A MMC scalar can only interact with gravitons through its kinetic energy but this cannot 
mediate any secular growth, in spite of the growing graviton field strength, because the scalar's kinetic 



energy redshifts to zero exponentially fast. By contrast, a fermion interacts with gravitons through 
its spin, in addition to its kinetic energy, and the spin-spin interaction remains effective even when 
the kinetic energy redshifts to zero [9]. The same thing seems to be true of a small mass [10]. 

The importance of spin in mediating interactions between inflationary gravitons and massless 
fermions suggests that there might be comparably strong effects on other particles with spin such as 
photons and gravitons. The graviton contribution to the one loop graviton self-energy has been worked 
out [11] but so far not used to quantum-correct the linearized Einstein equation. The purpose of this 
paper is to derive the one loop graviton contribution to the vacuum polarization on de Sitter back- 
ground. We will use it to solve the quantum-corrected Maxwell equation in a subsequent paper. (The 
flat space analog of this problem was carried out as a warmup exercise [12].) Our computation is done 
in dimensional regularization, fully renormalized with the necessary BPHZ (Bogoliubov-Parasiuk- 
Hepp-Zimmermann) counterterms [13], and reported in the noncovariant tensor basis whose efficacy 
has been demonstrated in a recent study [14]. We also use the Hartree approximation to argue that 
photons likely experience a secular suppression of their electric field strengths. 

This paper contains seven sections of which the first is this Introduction. Section 2 gives those of 
the Feynman rules of Maxwell -|- Einstein which are needed for our computation. The contribution 
from a single 4-point vertex is derived in section 3. Section 4 gives the much more complicated 
contribution from two 3-point vertices. Renormalization is accomplished in section 5. Although the 
use of our result to quantum- correct Maxwell's equation is deferred to a latter work, the Hartree 
approximation is employed in section 6 to argue that photons experience secular changes of the same 
strength but opposite sign as those of fermions [7, 9]. Our conclusions are given in section 7. 

2 Feynman Rules 

The purpose of this section is to present the formalism used to compute the one loop quantum 
gravity contribution to the vacuum polarization depicted in Fig. 1. We begin by describing the 
background geometry. Then we use the primitive Lagrangians to derive formal expressions for the 
first two diagrams of Fig. 1. The longest subsection discusses our conventions for gauge fixing and the 
resulting propagators. We next describe how the vacuum polarization can be represented in terms of 
two structure functions. The section closes by giving the counterterms needed for this computation. 

2.1 Our de Sitter Background 

We model primordial infiation as the cosmological patch of de Sitter space. The invariant element is, 

ds^ = a^ [—dr]'^ + dx ■ rfx] , (5) 

where a{ri) = —^ = e^^ is the scale factor and H is the Hubble parameter. Whereas the spatial 
coordinates x take their usual values, the conformal time r] runs from rj — )• — oo (the infinite past) to 
?7 — )■ 0~ (the infinite future). 

In representing functions such as propagators which depend upon two points, x'^ and x"^, we will 
make extensive use of the de Sitter length function, 

y{x]x') = a{r])a{r]')H'^ [\\x — x]]"^ — (\r] — r]'\ — i5Y~\ . (6) 

We also need the de Sitter breaking product of the scale factors a at x^ and a' at x"^, 

u = In(aa') . (7) 
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Figure 1: Graviton contributions to the one loop vacuum polarization. Photon propagators are wavy 
and graviton propagators are curly. 

Derivatives of y and u furnish a convenient basis for representing bi-vector functions of x^ and x'^ 
such as the vacuum polarization, 
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It turns out that either taking covariant derivatives of any of the five basis tensors (8), or contracting 
any two of them into one another, produces metrics and more basis tensors [15, 16]. 



2.2 Our Primitive Diagrams 

The total Lagrangian consists of the primitive contributions from general relativity and electromag- 
netism, plus the BPHZ counterterms necessary for this computation. 
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The primitive Lagrangians of general relativity and electromagnetism are. 
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The symbol G stands for Newton's constant, while A = {D — l)if^ is the cosmological constant. We 
employ a L>- dimensional, spacelike metric gf^,^, with inverse g'^'^ and determinant g = det(5'^,y). Our 
affine connection and Riemann tensor are, 
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Our Ricci tensor is R^u = R^'^p^ and the associated Ricci scalar is i? = g^^R^u- The electromagnetic 



field strength tensor and its first covariant derivative are. 
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We define the graviton field h^u{x) as the difference between the full metric and its de Sitter 
background value a?"ri^y, 

g^y{x) = a^{'q) Tj^y + Kh^y{x) = a^'g^,y{x) , (14) 

where h? = 16ttG is the loop counting parameter of quantum gravity. We follow the usual conventions 
whereby a comma denotes ordinary differentiation, the trace of the graviton field is h = rj^^h^j^p, and 



graviton indices are raised and lowered using the Minkowski metric, h^^ = rj^'^hpy and h^^ = rj^'^ri'^^hpa- 
Up to a surface term the gravitational Lagrangian can be written as, 
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From Fig. 1 one can see that we only need (15) for the graviton propagator. 

The only interactions we require descend from the second variational derivative of the electromag- 
netic action. 
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The necessary vertex functions are obtained by expanding the metric factors. 
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The tensor factors for the 3-point and 4-point vertices are, 
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Parenthesized indices are symmetrized and indices enclosed in square brackets are anti-symmetrized. 
If we call the graviton propagator ilap^-ysl^x] x') and the photon propagator i[pAo-](x; x'), we can 
give formal expressions for the first two diagrams of Fig. 1. The one constructed from a single 4-point 
vertex is, 
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The diagram constructed from two 3-point vertices is. 
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2.3 Our Propagators 

The quadratic part of the gravitational Lagrangian (15) is. 
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Before fixing the gauge and giving the graviton propagator we must digress to summarize the long 
and confusing debate between cosmologists and mathematical physicists concerning the de Sitter 
invariance of free gravitons [17, 18]. Although the propagator equation can be made de Sitter invariant 
by an appropriate choice of gauge, that does not guarantee the de Sitter invariance of the solution. The 
classic counter-example — which plays an important role in our solution for the graviton propagator 
— is the propagator iAA{x; x') of a massless, minimally coupled scalar, 

^A^f^g^''^JAA{x■,x')) = y/^ UiAAix; x') = i5^{x-x') . (23) 



Equation (23) is de Sitter invariant, but there is no de Sitter invariant solution for iA^(a;;x') [19]. 
This can be seen from the time dependence of the coincidence limit [20], 

iA^(x;x) = (Divergent Constant) + — ^ x ln(a) . (24) 

If one chooses the "E(3)" vacuum [21] to preserve the spatial homogeneity and isotropy of cosmology 
then the unique solution is [22], 



iAA^x; x') = iAc{{x; x' 
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where iAcf{x;x') is the (de Sitter invariant) propagator of a conformally coupled scalar, 

ttD—2 . r\ . / A\ — — 1 

tA^,{x-x') = -—j,T(--l)(-y . (26) 

(47r)T V2 J\yJ 

The massless, minimally coupled scalar is especially relevant to gravitons because Grishchuk 
showed that the physical components of h^y (transverse, traceless and purely spatial) obey the same 
equation [23]. Hence the graviton propagator must break de Sitter invariance as well. A cosmologist 
would also see this from the scale invariance of the tensor power spectrum [18]. Mathematical physi- 
cists for years disputed that conclusion because they found manifestly de Sitter invariant solutions to 
the propagator equation which results from adding de Sitter invariant gauge fixing functions to the 
quadratic Lagrangian (22) [24]. 

The discordant viewpoints have recently converged somewhat with the demonstration that there is 
an obstacle to adding invariant gauge fixing functions on any manifold, like de Sitter, which possesses 
a linearization instability [25] . That still leaves open the possibilities of either adding a noninvariant 
gauge fixing function or else enforcing a de Sitter invariant gauge condition as a strong operator 
equation. When either possibility is pursued, along with the requirement that the resulting propagator 
can be expressed as a superposition of plane wave mode functions [26, 27, 28], the result is a de Sitter 
breaking solution of precisely the form implied by the scale invariance of the tensor power spectrum 
[29, 30]. However, when a de Sitter invariant gauge condition is employed, in conjunction with analytic 
continuation techniques (either from Euclidean de Sitter space, or in the mass-squared of certain scalar 
propagators), the result is de Sitter invariant [31, 32], except for certain discrete choices of the gauge 
condition. 



It is well known that analytic continuation fails to recover power law infrared divergences [33, 34, 
16]. In this context the discrete problems, which have long been noted [35, 36], seem suspiciously like 
the special values at which a power law infrared divergence — which is always present — happens to 
become logarithmic and hence visible to an analytic continuation technique. We will therefore employ 
a de Sitter breaking graviton propagator, which seems to be the safer choice. Despite the continuing 
disagreements, it is important to note that very little difference remains between cosmologists and 
mathematical physicists. In particular, the various de Sitter breaking solutions for the graviton 
propagator all give the same result for the linearized Weyl-Weyl correlator [37, 38], and this agrees 
with the result from de Sitter invariant solutions [39] once some errors are corrected [40]. 

We fix the gauge by adding [26], 
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The resulting graviton propagator can be expressed as a sum of constant tensor factors multiplied by 
scalar propagators [26], 
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We have already seen the de Sitter breaking A-tjpe propagator (25). The 5-type and C-type propa- 
gators obey the equations. 
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Each of these propagators is de Sitter invariant and consists of iAcf{x;x') plus an infinite series of 
less singular terms which vanish in D = 4 dimensions. 
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Note that the -B-type and C-type propagators agree for D = A dimensions. The tensor factors are, 
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where tj^'^ = tj^^ + 5q5q is the spatial part of the Minkowski metric. 



The quadratic part of the electromagnetic action is, 
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Of course it is no more possible to add a de Sitter invariant gauge fixing function for electromagnetism 
than it is for gravity [25]. Unlike gravitons, photons show no physical breaking of de Sitter invariance, 
so the photon propagator is manifestly de Sitter invariant if one employs an exact de Sitter invariant 
gauge condition [41] . However, the de Sitter breaking of the graviton propagator implies that there is 
no point to keeping the photon propagator invariant. We have chosen instead to add the noncovariant 
gauge fixing function which is most closely related to the gravitational one (27) [42], 
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The associated photon propagator is [42], 
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2.4 Our Structure Functions 

The vacuum polarization is a bi-vector density i[^'n.'^]{x; x') which is transverse at each point. 
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Although it possesses 16 components in D = 4 dimensions, the combination of transversality (38), 
refiection invariance — [TI'^] (x; a;') = ['^11''] (x'; x) — and the coordinate symmetries of the vacuum 
relate these components so that they can be expressed in terms of a very few structure functions. For 
example, Poincare invariance implies the following simple form in the fiat space limit. 



tiTTiy 



W 



flat 



^lu pa _ 



ri'"'ri''ndpd'^n{x- 



-X 



(39) 



The one loop fermion and scalar contributions to Il{x — x') have been known for decades, and an 
explicit result for the one loop graviton contribution has recently been derived [12]. 

Because the graviton vacuum on de Sitter background does not respect full de Sitter invariance, 
but only spatial homogeneity and isotropy, it turns out that two structure functions are needed [14]. 
We could still choose to represent the transverse projection operators using covariant derivatives. 
However, a detailed examination of this form for the already derived vacuum polarization from SQED 
[3] reveals that it is cumbersome and that it obscures rather than simplifies the essential physics [14]. 
This seems to be because the conformal invariance of classical electromagnetism in D = 4 dimensions 
is a more powerful organizing principle that the background's isometrics. We have therefore chosen 
to employ the noncovariant form originally used to represent the SQED result. 
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We remind the reader that fj^^ = r]^^ + 5q5q is the purely spatial part of the Minkowski metric. 

By comparing (40) with (39) one sees that our structure function F{x; x') must agree with the 
flat space result iH(x — x') in the limit that H vanishes with the co-moving time t = ln[a)/H held 
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fixed. Hence the leading divergences are contained in F{x;x'). All terms in G{x;x') must contain 
at least one factor of H^, and they are correspondingly less divergent. Although our representation 
(40) is not de Sitter covariant, the two structure functions have very simple physical interpretations 
in terms of changes in the electric permittivity and the magnetic permeability [43]. However, there 
is a straightforward procedure for converting our results for F{x; x') and G{x; x') to the physically 
opaque, de Sitter covariant representation [44] if that is desired. 

2.5 Our Counterterms 

Deser and van Nieuwenhuizen showed that Einstein + Maxwell is not renormalizable at one loop 
order [45]. However, it is straightforward to absorb the divergences, order by order, with BPHZ 
counterterms [13]. We can then solve the quantum-corrected Maxwell equation (3) in the standard 
sense of effective field theory [46] . This has already been done for quantum gravitational corrections 
to electrodynamics on flat background [47, 12]. 

The vacuum polarization has two external photon lines so our computation requires counterterms 
with two vector potentials. The superficial degree of divergence is four at one loop order, which means 
there must be four derivatives acting either upon vector potentials or metrics. However, U{1) gauge 
invariance implies that at least two of the derivatives must act on vector potentials. If both of the 
remaining two derivatives act upon vector potentials we have the single counterterm which survives 
in flat space [12], 

C^D^F,,DpFp„g''Pg^Pg''''./^g . (41) 

Because our gauge fixing functions (27) and (36) reduce, for iJ — )■ at fixed co-moving time, to those 
employed in our previous fiat space computation [12], the divergent part of C4 must agree as well, 
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There are three invariant counterterms with two derivatives acting on two vector potentials and 
two acting on metrics, 

CiF^uFp^g^''g''^R^^+C2F^.F^„g^^^m''"^/^g + C^F^,Fp„R^^''p'^^^g. (43) 

However, all the curvatures are related in de Sitter background, 

(44) 

(45) 
(46) 

Our computation therefore determines only the combination C = (D — 1)-Dx Ci + (L' — 1) XC2 + 2XC3, 
and we can write the resulting counterterm as just CH'^F^^Gpag^^g^'^ \/—g- 

Because our gauge fixing functions (27) and (36) break de Sitter invariance, it is necessary to 
consider counterterms which are U{1) invariant but not generally coordinate invariant. Two properties 
of our gauge conditions restrict the number of possibilities: 

• They become Poincare invariant in the fiat space limit of iJ — )■ at fixed co-moving time; and 

• They preserve spatial homogeneity and isotropy, as well as the dilatation symmetry x^ — )■ kx^ . 
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The first property means that (41) is the only counterterm without exphcit factors of H. The second 
property imphes that the only extra counterterm we require has the form, d^CH'^FijEk^g'^'^g^^^—g. 
Our computation therefore requires only three counterterms, 

£bphz = ^CH^F,,Fu,g'^g^'^^g + CH^F^,Fp,g>'Pg'"'^^g 

+C^D^F^,DpF^„g^^g>^pg^'^^^g . (47) 

It remains to work out how the counterterms (47) affect the structure functions F{x; x') and 
G{x\ x') in our representation (40) of the vacuum polarization. The first step is taking the variation 
with respect to A^{x) and A^i^x'), 
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The next step is specializing to de Sitter, with g^^u = a^rj^i, and V^^^^ = aH{6j^6^ + 5^5° — rjp'^ri^^). This 
is very simple for the AC and C counterterms. 
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Of course (50) contributes directly to the structure function F{x; x'), and (49) contributes to G{x; x'). 
The Ci counterterm is complicated because of the way the tensor d'Alembertian acts on F„^, 
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The first term of (51) contributes only to the structure function F(x;x'), whereas the second and 
third terms contribute to both structure functions. After some tedious tensor algebra and reflection 
of derivatives on the delta function {8^5^ {x — x') = —daS^{x — x')), we find that the counterterms 
make the following contributions to the two structure functions. 
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The simplicity of this form is one more indication, added to the many already found [14], of the 
superiority of the noncovariant representation over the covariant one. 



I 


jjUpK.XaP'^S 


1 


1 k[A p]/i a/3 7<5 


2 


_ 1^k[A p]^ a(7 <5)/3 


3 


Q/3 7)[/i k][p A]{5 _|_ ^'y5^a)[fi K.][p X](p 


4 


rpi^ia I3)[\ p\(^ 5)p _|_ ^k{7^5)[A p](a /3)/i 


5 


K(a /3)(7 5)[A p]/i _|_ ^K{7^5)(a /3)[A p]p 


6 


j^k[X p]{a I3){y S)p _|_ r^^KlX p\(^ 5)(a I3)p 



Table 1: Parts of 4-point vertex function 

3 The 4-Point Contribution 

This section summarizes the derivation of the contribution to the vacuum polarization from the 4- 
point diagram in Fig. 1, by deriving the structure functions required for the desired representation of 
the vacuum polarization (40). This is achieved by first completing all naive index contractions of (20), 
followed by a substitution of the graviton propagator, which allows all final index contractions to be 
taken. Lastly, several contraction identities are introduced from which the 4-point structure functions 
can be deduced. The procedure outlined in this section will be similar to that used for deriving the 
3-point contribution and serves as a simple guide to extracting the desired structure functions. 

3.1 Naive Contractions 

In order to derive the two scalar structure functions from the 4-point contribution to the vacuum 
polarization (20) the first step will be to substitute (19) for the 4-point vertex function and apply the 
naive contractions. This will naturally engender many terms, and to make the process simple for the 
reader to follow we will break (19) into pieces and then define the 4-point contribution as a sum of 
those pieces 



p-nv 
■'-'■4pt 



(x;x') = 9j2K2a^-^^f/r"^ 



afi-yS 



1=1 



a(5^-y5 



{x;x)dx6 {x — x') 



(54) 



where Table 1 lists the various f/j:""' "'^ '^ . Once the vertex function terms have been inserted into 
the vacuum polarization the naive index contractions can be carried out term by term, the results of 
which are listed in Table 2. 
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I 


9Jzfi:2ai3-4f/M<^-Aa/375^" 


aP^-lS 


{x-x)dx5''{x-x')\ 


1 


Ir^^^d, [in'a^-H 


'' ,A^ ; 


5«5o|_i5-|,;,2^o-4, 


' ,A^ ; 


9^5^} 


2 


-Iv^^'d^ [iH'a^-S 


^'A,; 


5«5^| + l5-|,;,2^^-4, 


"^'A,5 


9^(5^} 


3 








-4. 

"/9 




4 


d, l^/t^a^-S 


liu \n\ 


d^6'']-d,[tK^a''-H 


A/1 A KI/ 


9a5^} 


5 


r^'^'^d^ iiK^a^-H 


XS^n ^ 


dx5'']-d^[iK'a''~H 


i.5^«^ 


9^5^} 


6 


d, ^m'^a^-H 


1.5 ^u ^ 


' d^5''^-d^[iK''a^-H 


^^A\ 


"9a5^} 



Table 2: Terms of 4-point contribution after naive index contractions. 

3.2 Substitution of Graviton Propagator 

To complete the index contractions the full graviton propagator must be inserted. This will again 
create many more terms so it is useful to break the graviton propagator into three pieces and consider 
each part separately. Upon consideration of the graviton propagator (28) we see that if each of the 
three types of scalar propagators are set equal, to say B{y), then the tensor components combine to 
give the conformal graviton propagator tensor component 



rpA 

l^V''' pa 


+ 


l^v-'- pa 


+ 


rpC 

l^y-'' pa 


= 


rpci 

py-'- pa 



'^Vp{pVa)y - 



{D-2) 



VpvVpa 



(55) 



By adding and subtracting B{y) from each of the scalar propagators in (28) the graviton propagator 
can be rewritten as 



pv'—^pa 



cf 
fJ-i^-^ pa 



t: 



X B 



rpA 

t^f-'- pa 



X HAa -B) + 



rpC 

f^^-^ pa 



x{C-B) . 



(56) 



Writing the graviton propagator in this way cancels the conformal parts of the scalar propagators 
in the second two terms. This property is not useful for the 4-point contribution, but it will be helpful 
for renormalizing the 3-point contribution and we use it again here for consistency. We will refer to 
the three components of (56) separately as the conformal part, the A-type part and and C-type part 
as associated with the tensor factor of each piece. 

Before making substitutions for the graviton propagator we need to know how each of the tensor 
components contract in the various ways appearing in Table 2. The relevant contractions are listed 
in Table 3, where each element represents the tensor factor in the top row contracted with the 
combination of fiat space metrics listed in the left hand column. 
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rpci 



rpA 

1^1^^ pa 



rpC 

iJ'V-l- pa 



^pu^pa 



^pp^ua 



^pu 



pp 



4D 



iD-2) 

D{D^-D-4) 
iD-2) 



4(g-l) 
(D-3) 

{D^-AD'^+D+2) 
(D-3) 



4 4 

'(Z)_2)^P'^ ~(D-3)^P'^ 



(D^-D-4) (D^-3D-2) - 

(D-2) ^^'^ (D-3) ^^'^ 



(Z)_2)(D-3) 



, (D^-5D+8) 
'(D-2)(D-3) 



4 .^0^0 , 4 



(D-2)"p"<T + {D-2)(D-i)'^P" 



^(D-2)"i^"f7 ^ (D-2)(D-3)'''^'^ 



Table 3: Various graviton tensor factor contractions. 



As stated above we dissect the graviton propagator substitution into three parts. The substitution 
and following index contractions for the conformal part of the graviton propagator are listed in Table 
4, the terms resulting from the A-type part are listed in Table 5, and the terms from the C-type part 
in Table 6. Note that in these tables we suppress the {x — x') factor on the delta functions for brevity. 
Upon further inspection of the conformal part of the graviton propagator in Table 4 we see that all 
six sets of terms have the same tensor structure. Combining all of the terms we find the simplified 
expression 



PY[v 



n 



4pt 



I iXy > lAj 



(1)3 - 9D2 ^ lOD + 16) 



cf 



-d'' [iK'a''-^ 



A{D - 2) 

B{0)d^'6^{x - x')] } 



{r]^''d, [iK^a^~^B{0)d''6^{x - x')] 



(57) 



Equation (57), and Tables 5 and 6 make up the 4-point contribution to the vacuum polarization. 
Next we will transform these results into the desired manifestly transverse form. 



3.3 Finding the 4-Point Structure Functions 

Recall that it is our goal to write the result for the vacuum polarization in the form of (40), where 
two transverse projection operators are acting on two scalar structure functions F{x; x') and G{x; x'). 
Using the Tables of section 3.2 we can now derive these structure functions. 

The easiest way we found to extract the structure functions was to exploit the known transversality 
of the vacuum polarization and isolate the structure functions via two contractions, resulting in two 
equations for the two structure functions. The first contraction is made with 5° 5°, applied to equation 
(40) this gives the identity 



OttO 



n^ 



(x: X 



V^F{x]x') , 



(5^ 



which provides a simple equation for F . To find G the second contraction is taken with Sl^6l{j ^ i), 
applied to (40) we find the second identity 



'W 



(x; x') = -d''d^ [F{x- x') + G{x- x')] . 



(59) 
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2^D-4r/-/^'^«:Aa/37<5 



dA iK'a^'-^U'i 



rpci 



B{0)dx6^{x-x' 



2(D-2) 



^ {r]'"'d^ [iK^a^~^Bd^6^] - d^ [iK'^a^~^B^^'8^] ] 



D{D'^-D-4) 
4{D-2) 



{r7^^(9« [iK,'^a'^"^Bd^6'^] - d" [iK,^a'^-^Bd''6^] } 



(D-2) 



{r^^-S, [ 



iK^a^ "' 



Bd'^S^] - d" [iK^a^-^Bd^d^] } 



^ {r]>'^d^ [iK^a^-^Bd^5^] - d" [in^a^-^Bd^'d^] ] 



(D-2) 



(Z)2-Z)„4) 
(D-2) 



{r]^'''^^ [iK'^a^-~*Bd^5^] - d" [iK^a^-'^Bd'^d^] } 



(D-2) 



{r]^"'^^ [iK.^a^-~^Bd'^5^] - d'' [iK^a^-^Bd^d^] } 



Table 4: 4-point contributions coming from the conformal part of the graviton propagator. 



2 D-4r/-M'^«^"/37<5 , 



9«<^2/t^a^-^f/; 



rpA 



(i^A{x-x)-B{Q)) dx6^{x-x' 



2(D-3) 



{r]^"'d^ [iK^a^~\iAA - B)d''5^] - d" [iK^a^^^iAA - B)^^'5^] } 



(D3-4D^+D+2) 
4(D-3) 



{r]^'''^^ [iK^a^-\iAA - B)d''6^] - 9^ [i^'^a^-^iAA - B)d^5^] } 



^ {V''9. [ 



(D-3) 



tK^a^ ^[t 



Aa - B)d^5^] - d" [iK^a^-\iAA - B)d^5^] 



+f]^"'d^ [if^a^-\iAA - B)d^5^] - d"" [iK'^a^-\iAA - B)^^'5^] } 



(D-3) 



{r7^^9« [iK^a^^\iAA - B)d''5^] - d" [ 



IK a [I 



Aa - 5)9^^5^] } 



'-'^V^ {^-^« [ 



iK?a^ '^[i 



Aa - B)d^6'^] - d" [iK-a-^-[i 



-2^D-4/ 



A^ - B)^^'5'^] ] 



(D^-3D-2) 
(D-3) 



[r^'d. [ 



iK^a^ ^[i 



Aa - B)d^6^] - d" [iK^a^-\iAA - B)d>'6^] } 



Table 5: 4-point contributions coming from the A-type part of the graviton propagator. 



13 



d^{iK^a^-^Ur^"^'^' I a^T^s {C{0)-B{0))d,6^{x-x') 



-2)(D-3) W ^'^ [ 



{D-2)(D-3) 



iK^a^ ^ 



{C - B)d^6^] - d'^ [ 



iK^a^ ^ 



(C - 3)8^^5^] ] 



2(D-2)(D-3) V' ^'^ L^'^ " 



(C - B)d^5^] - d" [iK^a^-\C - B)d^'6^] } 



(D-2) 



f^f^ogu \^i^2^D-i^^ _ 5)^0 ^Dj ^ j^^d" [iK^a^-\C - B)d''6^] 



_^M^50 [^^2^D-4(^ _ B)d°5^] - j^v^'^d^ [iK^a^-^ 



{C - B)d''6^] 



7]^^dP [iK^a^-\C - B)d>'5^] + T^d" [iK^a^~\C - B)^^'5^] 



.^mo^^o^^ [m2a^-^(C - B)d^5^] - j^^v'^^'d^ [iK^a^-\C - B)d^5^] 



{D-2)(d-3) 



{{D - 3)f]^''d^ [ 



2^D-4 



IK a 



{C - B)d^5^] + f]^"'^^ [iK^a^-\C - B)d''5^] 



-{D - 3)r7^V°'9« [ 



in'^a^ ^ 



{C - B)d^6^] -{D- 3)r]''^d^ [ 



m'^a^ ^ 



{C - B)d^b^\ 



-{D - 3)?7^°(9'^ [ 



iK^a^ ^ 



{C - 5)^0(5^] - d'' [iK^a^'\C - B)d''6^] } 



(D-2)(D-3) 



{-(D- 3)^-9° [ 



2^D-4 



tK a 



Dl au \-^2^D~i 



{C - B)dH^] - d" [ 



IK a 



\C - B)^^'6^] 



{D - 3fr]''°d^ [it^a^-\C - B)d^'5^] + T]^"'d^ [ 



iK^a^ ^ 



(C - B)d^6^] } 



{D-2)(D~3) 



[n^^'d^ [m2a^-4(C - B)d^5'^] - d" [iK''a^~\C - B)d^'5^] 



-{D - 3) VV^-^^ [iK^a^~\C - B)d''6^] + {D - 3f^]^'^^'' [iK^a^-\C - B)d^6^] } 



Table 6: 4-point contributions coming from the C-type part of the graviton propagator. 
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Knowing F it is now trivial to find G. 

It is true that there are other pairs of contractions that would work equally as well. For example, 
contracting with 5^51 and fj^^ leads to the identities 






{X] X 

(x; x' 



-d"'d'F{x; x') , 

-{D - l)9°9"'F(x; x') -{D- 2)V=' [F(x; x') + G(x; x')] 



(60) 
(61) 



These identities provide a useful check on the structure functions derived from identities (58) and 
(59). There are likely more contractions that would provide similar sets of equations, but these two 
sets were sufficient for our purposes. 

It is quite a tedious task to show how applying these contractions plays out for the entire 4-point 
contribution, and the procedure is extremely repetitive. For the reader's sanity and our own we will 
work out one example and assume that the procedure for the rest of the 4-point contribution can be 
easily deduced. We will demonstrate how to find F and G from the conformal part of the graviton 
propagator as given in (57). 

Contracting 5^5^ with (57) we have 



OttO 



n 



4pt 



Jcf 



X] X 



{D^ - 9D 



4(Z,-^^{-^-[-'''°"'^^"^"' 






{D^ - 9D^ + lOD + 16) 
A{D - 2) 
{D^ - 9D^ + lOD + 16) 
4(D-2) 



-d'[tK'a''-^Bd%''{x-x')]} 
[{d-d' + d'^d"^) X [m2a^-^55^(x - x')] } 

V^[tK^a''~*B{0)6''{x-x')] , (62) 



where in going from the first line to the second we used the delta function to make the change d — ?> —d'. 
It is then trivial to pull the inner derivative outside the curly brackets since all prefactors are only 
functions of x. Also the B-type propagator can be evaluated at y = since it is being evaluated at 
coincidence {x = x'). Upon comparison with identity (58) we find the first structure function 



Vfe.-)='^^'-^5^^°^ + ^°'[.Va°-- 



4{D-2) 

To find the second structure function we contract Sl^6i{j ^ i) with (57) 

{D^ - 9D^ + lOD + 16) 



5(0)(5^(x - x')] 



(63) 



'n^ 



[x: X 



cf 



4(D - 2) 
{D^ - 9D^ + lOD + 16) 



jo - d^ [if^a^-*Bd'6^{x - x')] I 
d"d'[i^a''-^B5''{x-x')] . 



A{D - 2) 

Invoking identity (59) and making the proper substitution for F4 ^ it is easy to see 

GA,cf{x;x') = . 



(64) 



(65) 



This concludes the example case for finding F and G. The same procedure can be applied to all 
of the terms in Tables 5 and 6 to find the rest of the structure functions. Combining all F and G 
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contributions from the three parts of the graviton propagator produces the full result for the 4-point 
structure functions 

F,ix-x') = .^a^-^{^^^^AA0)-^-^Bi0)-^j^^C{0)}^6^{x-^ (66) 

G,{x;x') = -n'a''-'{^-^^j^^^tAAiO) - iD-3)BiO) -2^^^CiO)y5''{x-x') (67) 

We can now make substitutions for the propagators. The B- and C-type propagators are a finite 
constant at coincidence in i5 = 4, but the A-type propagator is divergent. Therefore, it is useful to 
break the structure functions into their finite and divergent pieces 

i^4,div(x;a:') = ^— '- ^ ^^,^ \, ^vrcot (-DJ ^^^(x - a;') , (68) 

F4,iimteix; x') = -—- <j- ln(a) > i6'^{x - x') , (69) 

G,M^;x') = (^_3^ S;0^Vif(|r^"^°H2^)^^ (^-^)' (^°) 

/.2^2 r I ^ 

G4,fimte(a;; x') = — — - <^ - + ln(a) > i5'^{x - x') . (71) 

This concludes our derivation of the 4-point contributions to the structure functions, where we note 
that only the terms proportional to ln(a) in (68-71) cannot be absorbed into the counter terms. We 
will now derive the 3-point contribution. 

4 The 3-Point Contribution 

This section will cover the derivation of the contribution to the vacuum polarization attributed to the 
3-point diagram in Fig. 1. The procedure for deriving the scalar structure functions from the diagram 
is very similar to the one used in section 3. First the naive index contractions are completed in pieces 
by dividing the 3-point vertex function appropriately. Then once the proper substitutions have been 
made for the graviton and photon propagators all remaining index contractions can be completed. 
Finally the 3-point contributions to F and G will be presented in several parts, broken up according 
to which pieces of the photon and graviton propagators the structure function originated from. 

4.1 Naive contractions 

Following the same organizing procedure as in the previous section, we complete the naive index 
contractions first by breaking the 3-point vertex function (18) into pieces, shown in Table 7, and 
rewriting the 3-point contribution as a sum of these terms 



■'■■'■Spt 



1=1 



(x: x') 



2 
,/Z)-4 ' 



xtna 

j=i 
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J2vr"^'''d,d'^z[,A^]{x;x')]. (72) 



T//ipKAa/3 



^al3^K[X^p]fi 



4^a)[M^K][p„A]{/3 



Table 7: Pieces of 3-point vertex function. 
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J 


-d^d'JK\aay-%'''^^''^t 


a/S'^-yS ( 


x;a:')V7''^"'9,9;^ 


pA, 


(x;x')} 


1 


1 


-K'd^d's{{aa')''-H 


'' A' s 


9[^9'[['^i 


K]^e]] 


} 


1 


2 


K'd^d's{{aa')''-H 


fi ^A<PW 


Qll^^Q'^, p 


lA^]" 


-(aa')^-S[ 


3 ;3A-[^ 


/ ^[[M^'e], 


"«iia; 


} 


2 


1 


K'd^d's{iaaY-H 


A[^^5 ^ 


^A^'ll'^^f 


^lA^ll 


-(aa')^-'^ 


p[^^5 


sld^^d'^'^i 


pA^l" 


} 


2 


2 


-K'd,d',{{aa' 
+ {aa')^-H 


p[H^O 


"a[p 


gn]gn 


5a 


s;^ 

^pA 




]]j-(aa') 




Q. 

9]], 


«]A 


p 


A^ll 

} 


- 



Table 8: 3-point results after naive index contractions. 

(72) contains products of the parts of the two 3-point vertex functions; the result for completing the 
naive index contractions for all said products are shown in Table 8. To simplify the results in this 
table a short hand notation for antisymmetrization has been adopted, where both square brackets 
and double square brackets imply antisymmetrization. The index structure for these terms can be 
confusing and it should be noted that the antisymmetrization only applies to the immediate indices 
with in the brackets. Here is a worked out example 



p[p^cr[[p 



QKlgn 



Pfi^cru 



d^d"^ - i 



pK A (TV 



'A 



Q^Q'e 



PK \(7 



A' 



d^d'" - i 



pp^crO 



To complete the index contractions, substitutions for the graviton and photon propagator must be 
made. 



4.2 Graviton and Photon Propagator Substitutions 

To make substitutions for the graviton and photon propagators and completing the index contractions 
as clear as possible we will be using the new form of the graviton propagator (56) and also break the 
photon propagator up in a similar manner. To modify the orginal photon propagator (37) we can 
again add and subtract B{x;x') from each term. Rearranging we find 



pA, 



I tj-/ • tXj 



')=r^,^aa'B{y) + 5l5l[B{y)-C{y)] . 



(74) 
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Since we will be considering different parts of the propagators individually for the rest of this section 
it is necessary to take a moment to explain our notation. We will refer to the first, second, and third 
term of (56) with the subscripts B, A, and C respectively. Likewise, for (74) we will refer to the first 
and second terms with the subscripts B and C respectively. 

In the 3-point contribution it is always a product of the graviton and photon propagator that 
appear, so there are in general six combinations of terms that can arise BB, AB, CB, BC, AC, and 
CC, where the first letter stands for the part of the graviton propagator being considered and the 
second letter stands for the part of the photon propagator. However, it can be shown that the last 
two combinations do not need to be calculated since they are made of the product of two differences 
of scalar propagators. For these cases the conformal parts of the scalar propagators will cancel and 
what remains are only the infinite sums in each propagator, but these terms go like ~ y^. Their 
degree of divergence is such that we can take D = 4 immediately, and in this case all of the infinite 
series vanish. It is true that the A-type scalar propagator has an extra term going like ~ y^~^l'^\ 
however, in the D = 4 limit the infinite series multiplying this extra term will cause it to vanish too. 
So, in the end there are really only four combinations that need to be calculated. 

One other notational comment is needed before the calculation can continue. All four possible 
combinations of the parts of the propagators will result in terms of the same form, namely 

#a„9;,{(aa')''"'^A,'9^9Kaa'?A,)} , (75) 

where il^x.y is actually a difference of propagators for the cases x = A,C or y = C. Since all of 
the terms will look almost identical, with only the propagators and indices changing, we can vastly 
simplify reporting our results by adopting a notation of only writing out the derivatives and their 
associated indices. These are the only parts needed because it is the indices that will make up 
the primary difference in each term, and the specific propagator combination can be denoted in the 
subscript of the vacuum polarization. In this notation (75) would take the simple form dad'p{d^d'^) ■ 

We are now ready to dive into the calculation. First we will perform the substitutions for the B 
part of the photon propagator with the B part of the graviton propagator. All of the terms for this 
portion will take the form 

i^dc,d'p{{aa')^-^Bd^d's{aa'B)] . (76) 

The full result for this set of propagator pieces is 



■'■■'■Spt 



(L)2 _ 4£) + 2) 



[x- x') = ^— - — TT^i^H -r7^^9„9;j(9"9'^) + d'd'^id^'d'^) - d ■ d'id'^d" 

B,B \^^ ^j 



+do,d'''{d"d"' 
+K^l-2ri^"'d ■ d\d ■d')-d- d'id'd'^') + 2d''d'>'{d ■ d') + 9«9'^(9'^9'") 

-7]^''d^%{d^d"') + d''%{d^d'^'\ . (77) 

Next we compute the vacuum polarization from the B part of the graviton propagator and the C 
part of the photon propagator. These terms will all take the generic form 

^d^d'f^ {{aa'f-^Bd^d's{aa'{B - C))} . (78) 



The full result is 



-^3pt 



(x: X 



B,c 



^2 r MO^i^O^^^/^^^a^//?) _ r^'^OgOg^^^gt.gffi^ ^ d^d'^d^d'") 



P-2) 
-r]''''dad'\d''d"'] 

+d ■ d'id^d'") + r]''^d°d'>'{d ■ d') + 7]^'^d^d'p{d^d'^') - r7'^°9«9'^(9"9'°) 
+77'^°9 ■ 9'(9''9'°) + r7^°9 ■ 9'(9°9"^) - 9«9:,(9"9"^) - 9°9'^(9°9"^)} . (79) 

When we consider the A part of the graviton propagator and the B part of the photon propagator 
the terms all have the form 



i^d^%{{ 



aa 



l\D-Ai 



iAA-B)d^d'siaa'B)} 



30) 



and the full result is 



^3pt 



x; X 



A,B 



{D 1){D ^l^2(^.g^g^^^gag>f3^ ^ d''d'^{d''d'^) - d ■ d'{d^d'' 



P-3) 
+K^-ri>''d ■ d'{d ■d')-d- d'{d''d'^) + d^d^^id ■ d') + d^d'^id^'d 



-7] 



^^""8 ■ B'id ■ B') - T]^"'dJ'^{d^d"') + d'^B'^'id ■ B') + d'B'^ip^B'''] 



K^ld^d'^'id'^d"') - ri^'dj'^id'^d"^) - d ■ d'id^'d") + d'B'^id^'d'^) 



{D-3) 

+d ■ d'{B>'B"') - B^d'^'iB'^B") + ri>"'BJ'^{B''B"^) - d'B'^iB^'B'^) 

+r]'"'dJ'^{d''B'^) + B ■ B'id'^B"') - B'^B'pid^B'^) - dJ'f'id''B"') 
-B^d'piB^'d'^) + f]f'''Bad'f^{B''d'^ - BJ'''{B''d") + B ■ B'iB^d'" 



(D-3) 



The last case is the C part of the graviton propagator and the B part of the photon propagator. 
These terms take the form 



#9^9^ {{aa')^-\C - B)d^d's{aa'B)} , 



^2) 
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and the full result is 



■'■■'■spt 



[x: X 



2?^ — ^K^\-v''^do^d"'{d''d"') + r]'Y''dMd''d'^) + d^d'^d^d"") 






(iJ — 2) < 
2 



iD-2)iD-3) 



K^-r]^"'d^d'{d''d"^) + d^d'Jd^'d"^) - d ■ d'id^d'") + d^d'^id^d'' 






iD-2) 



-r]^'%''^do,B'^{d''B'^) - f]^''da,d'\d''d"^) - r]^'^B ■ B'id^d"") 

_^o^/0(^M5'-)| . (83) 



4.3 Finding the 3-Point Structure Functions 

Recalling the procedure outlined in section 3.3 we will now find the 3-point contributions to F and 
G in much the same way. The process will be a little more labor intensive since we no longer have 
a delta function on each term to assist in pulling out internal derivatives; there is also the added 
complication of having two internal derivatives instead of one. These changes will be accounted for 
as follows: 

First we can still use the same identities (58) and (59) as our set of equations for finding F^ and 
G4. Next we notice that in the absence of a delta function we can no longer simply change d — )■ —d', 
instead we have to carefully account for feed down terms that will arise from extracting derivatives. 
All of the 3-point terms take the form of (75), thus there are only four possible combinations of 
internal derivatives we will encounter, they are did',, dod'^, (9j(9q, and BoOq. Extracting these sets of 
derivatives will always result in the same feed down terms regardless of the propagators involved. A 
shght modification is needed for the A-type propagator, but in general extracting these derivatives 
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results in the following identities 



(aafSA^dod-iaa'iA. 



dod'i \{aa'f-^I^[iAjAy]\ + d'AHiaa'f^'^a (l[iAjAy] 



-{D - ?>)I^[iAjAy\ + /^[iA^iAj 



j/j ; f ' 



lA) 



{aa')^~HA^di%{aa'iAy 



dA 



[aa 



nD-4T2\ 



r[iAjAy] + dA H{aay--^a I[iAjA 



iD-3)P[tA,,iA''] + P[tKiAy] 



^5) 



{aa'f-HA^did'^{aa'iAy) 



d.d'^ 



aay-^r[iAjA, 



-2r]i^ 



H'{aaY~ I[i^xiA^ 



m) 



{aa')^-HA^dodo{aa'iA 



-2{D - 3)H\aa'f-^I[iAjAy] + dod'o {aaY-^r[iAjAy] 

+ {doa' + d'^a) x |i7(aa')^~^ I[iAjAy] - {D - 3)P[iAJA, 
1 



+ -{aaY-'V'r[i A J Ay] + H'{aay-\iAJAy + (2 - y)I[iAjA, 



-{D - l)I'[iA,,tAy] -{D- 3yi'[iA^iAy] , 



^7) 



where / represents and an indefinite integral with respect to y, and a dot over the propagators 
represents a derivative with respect to y^ . The A- type propagator is unique in that it is not just a 
function of y, but also contains a de Sitter breaking piece. It can be rewritten in the form 

iAA{x; x') = A{y) + ku , (88) 

where k = (^ \d/2 t(d~I2) ' ^"^ ^^^'^ form we can see that identities (84-87) will miss the feed down terms 
that arise when one of the internal derivatives act on the second term in (88). To account for this the 
following terms need to be inserted in each of the above identities when iA^ = {iA^ — B) 



\aa 



>\D-A 



kudodiiaa'B) -^ -d[ [H. 



a[aa 



r\D-3 



kB] 



^9) 



(aa 



/\D-4 



kudid'oiaa'B) -^ -di [Ha'{aaY~^kB] 



(90) 



aa')^"^kudid'Aaa' B) ^0 



(91) 



j\D-ii 



[aa:r ^kudod'^^aa' B) ^ -{d^a' + %a) x [H{aa')^'^kB] 

+2H\D - A){aa'f-^kB 



(92) 



We are now ready to apply the contractions in identities (58) and (59) and then use the substi- 
tutions (84-87), making use of (89-92) where appropriate, to find F^ and G3. Again, we will work 
through an example for the reader and then state the final result for all of the terms. 



^One might worry about local delta function contributions from terms of the form dodgB, but actually these contri- 
butions are fully accounted for in delta functions arising from feed down terms of the form OoOq {aa')^~^ PliAxiAy] 
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We consider all terms with the coefficient - — tjjZ2)^'^ ^^ i^''')^ since these form the smallest set of 
transverse terms. Performing the ffist contraction with (5°(5[^ results in 

_ {D^ -AD + 2) 
B,B ~ (D-2) 



OttO 



n' 



3pt,ex 



K^d^d'^ [{aaY-^Bd''d'%aa'B)] 



+d''d'p [{aa')^-^Bd''d'^{aa'B)\ -d-d' [{aa')^-^Bd''d"'{aa' B 
+dad'^ [{aaY-^Bd''d'\aa'B)]\ 
{D^ -AD + 2) 



{D-2) 
-V' [iaa'f-''Bdod'o{aa'B 



K^ d,d'^ [{aaY-^Bd,d'j{aa'B 



(93) 



where we have taken the 3+1 decomposition in going from the ffist line to the second. Now we can 
make the appropriate substitutions using identities (84-87) and compare with (58) to find 

i^3,cx(x;a:') = ^^'^^^_^^^ ^^ /.^{-2(D - 3)H\aaY~'l[BB] 



+ {doa' + d^a) x 
+H^(aaY-^ B 



aa')'^-H[BB] 
H{aa')^-'' {l[BB] - {D - 3)P[BB] 

2 I /A ^ \ T \ II>'^^ f T^ i\7'2rD2l { T^ o\2 7'2 r 



(4 - y)I[B'] -{D- l)r[B'] -{D- 3yr[BB] 



+ -[aa 



>\D-3t-,2 r3rD2i 



To find the companion G^^cx we contract Sli_6i{i ^ j) with the same set of terms in (77) 

JTTJ 
llq 



(94) 



3pt,ex 



B,B 



{D' AD + 2; ^2^ ^ ^ ^^.^, [{aa'f-'Bd'd'^aa'B 



{D-2) 



/\D-Auaa 



Bd'^d'Haa'B 



K^i-djd'^ [{aa')^-^B^^^'^{aa'B)\ + 8,8'^, [{aa')^-^B8,8'k{aa'B)\ 



-8 ■ 8' [{aa')^~^B8'8''{aa'B)\ + 8^8" [{aa: 

{D^ -AD + 2) 
- {D-2) 

+8o8'o [{aaY-*B8i8'^{aa'B)] + V^ [{aaY-*B8,8'j{aa'B) 

-8o8i [{aaY~^B8o8'j{aa'B)] + 8k8l [{aaY-^B8k8'^{aa'B)] , (95) 

where again we have taken the 3+1 decomposition in going from the ffist line to the second. Applying 
identities (84-87), (95) can be rewritten in the form 



«rrJ 



'n 



3pt,ex 



(^' - ^^ + \'8:8A{V' - 8,8', - 8',' - 8o')x 



,aa 



-d'o 



B,B {D - 2) 

'^""-H^IBB] 



-8o 



H{aa')^-'a' [I[BB] - {D - 3)l\BB] + l'[B] 
H{aa')^-^a (l[BB] - {D - 3)I^[BB] + I^[B] 



+AH\aa: 



l\D-2 



I[B' 



(96) 
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Substituting F^^cx into (59) we find 



^Sjexl^'C) -^ ) 



{D^ -AD + 2) 



iD-2) 
{D - 1)P[B^] - 2{D - 3)I[BB] - yI[B^' 



k" {H\aa')''-^ (b^ - {D - ?,fP[BB] 



1 



,/nD-3v72 r3rD2i 



+ -[ac^)—'V'r\B'' 



+ ((9oa' + d'^a) X 
+5' 



H{aa')^-^ U[BB] - {D - 3)P[BB] 
Hiaa'f^'^a (l[BB] - {D - 3)I^[BB] + I^[B^] 
Hiaa'f-'^a' (l[BB] - {D - 3)I^[BB] + I^[B^ 



+ {do + dl 



if 



[aa 



/\D-3t2 



r[BB] 



(97) 



This concludes the 3-point example; it should be clear how to proceed with the rest of the 3-point 
contribution. 

We now present the results for the scalar structure functions organized according to propagator 
combination. From (77) we find 



Fbb{x;x' 



{D^-D- 4) 
P-2) 



KHHHaa')''-^ 



2 [621 



B' -{D- 3yr[BB] -{D- l)r[B 



-2{D-3)I[BB] 
+H {doa' + %a) X 



[aci 



/nD-3v72 r3rD2 



[aci 



i\D-3 



I[BB]-{D-3)P[BB] 



[aa 



i\D-3t2 



F[BB] 



2tt2 



k'H 



{D-2) 
Gbb{x;x') = 



[aa: 



r\D-2 



[4D{2D -5)-{D^ -D- A)y] I[B^' 



{D-2) 
-2{D-3)I[BB] -yI[B'^ 



D{D - 3)H'^{aa' 
1 



/\D~2 



B^-{D- 3yP[BB] -{D- l)P[B^] 
D{D - ?>){aa')^-''W^f[B^] 



aa: 



-H {doa' + Sgo) X 

+ {D'^ -AD + 2)I[BB] 

-{D^ -4D + 2)Hdo 
-{D^-AD + 2)Hd'Q a'{aa' 



>\D-3 I t2iu2 



r[B^] - {D' -AD + 2){D - 3)F[BB] 



a(aa: 



f\D-3 



f\D-3 



^^L-^I'[B'] + {D- 3)I'[BB] - I[BB]^ 
^ P[B^] + {D-3)P[BB]-I[BB] 



+ {D^-AD + 2){do + dl 



if 



[aa 



{D-2) 



(9^ 



(99) 
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From (79) we find 

Fbc{x;x') 



(^{2P^-/^-4)/7 



2^aa')^-2 



I[B{B-C)] 



+ {3D - 8){aaY-'V'r[B{B -C)]\ , 



(100) 



Gbc{x'i x') = 2/t' 



XD-3) 



iD-2) 



DH\aa'f-'I[B{B - C)] 



[act: 



/\D-3t2 



I'[B{B - C)] 



:ior 



From (81) we have 

Fab{x;x) = K^lH^iD - l){aa)^- 



B{iAA -B)-{D- 3yr[{iAA - B)B] 



-{D - 1)P[{iKa - B)B] - 2{D - 3)I[{iAa - B)B] + (4 - y)I[{iKA - B)B] 
-r\{D - l){aa')''~'VH\iAA - B)B] 



+ (V^ + {D- l)d,d',) X 
+ {D - l)H{doa + 9oa) x 



,aa 



i\D-?,t2 



P[{iAa-B)B] 



,aa 



/\D-3 



IKiAa - B)B] -{D- 3)I\iAa - B)B] 



+ {D - l)kK^{^2{D - A)H\aa'f-^B - H{d^a' + 9» [{aa'f-^B] } , 

Gab{x- x') = -K^{H\aa'f-\D - 1) [^(^A^ - B) - {D - Sfl^HtAA - B)B] 

-{D - l)I^[{iAA - B)B] - 2{D - 3)/[(iA^ - B)B] 
(D-A) 



(102) 



-4 



(D-3) 



+ y I[itAA - B)B] 



+ 



{D^ -AD + 1] 



2 , Q'2^ 



v' + {D- 2){d', + d',') + {D- i)dod'. 



aaT I [i^^A-B)B] 



]^{D -l){aa')''~'VH\iAA- B)B] 



+H{D - 2)% 



a'iaa')''-'' ' ^' 



l'[{tAA-B)B]+I[{tAA-B)B] 



iD-3)P[itAA-B)i3] 



+H{D - 2)9o 



a[aa 



nD~3 I 7-2 



I'iitAA- B)B]+ I[{tAA- B)B] 



-{D-3)l\tAA-B)i3] 
-H{D-l){doa' + d'Qa) x 



,aa 



/nD-3 



mAA - B)B] -{D- 3)P[{iAa - B)B] 



-kK^h{D - A){D - l)H\aa')^-^B - {D - 2){doa + d'^a) x [(aa')^"^5] 



-H{D - l)d', [ 



I \Ji /\D-3 



a (aa 



B] - H{D - l)9o [ 



a[aa 



r\D-3 



B] 



(103) 



24 



And lastly, from (83) we have 



Fcb{x;x') = 2 



P-3) 
(D-2) 



K^H^iaa')''-' 



B{C -B)-{D- 'iYl\C - B)B 



~{D - l)P[{C - B)B] - 2{D - 3)/[(C - B)B] 



+ (A-y)I[iC-B)B] 
+H{doa' + 9oa) x 



-{aa 
2 



/\Z)~3v72r3 



V'r[{C-B)B] 



aa r " I I[{C - B)B] - {D - 3)r[{C - B)B] 



[aa 



/\D-3 
I\D-3t2 



r[iC-B)B] 



(104) 



Gcnix; X 



2k^ 



H^iaaY-' 



-{D - 3)B{C -B) + {D- 3yr[{C - B)B] 



{D-2 

+ {D - 3)(D - l)P[{C - B)B] + 2{D - d,fl[{C - B)B] 
{D-2){D-A) 



- 4 



P-3) 



iD-3)y I[iC-B)B] 



{D-A){D-2j^. ^^^ ^^^^, ^, ^„ 



P-3) 



-V^ + 2{D - 3)9o9^ - dl - d'. 



X 



aaT^'I [{C-B)B] 



1 



--{D - 3){aa')^~^V^I%C - B)B] 
+H{D - 3){doa' + d'^a) x \{aa' 



f\DS 



{D - 3)P[{C - B)B] - I[{C - B)B] 



-Hd. 



aiaa' 



i\D-3 



-I[{C-B)B] 



-Hd' 



aUaa')""-^ 



{D - 3)P[{C - B)B] - P[{C - B)B] 



{D - 3)J2[(C - B)B] - I\C - B)B] 



-I[{C - B)B] 



(105) 



Because of the different combinations of propagators appearing in each set of F and G it was not 
useful or enlightening to combine them here. Instead, once they are renormalized in the next section 
they are combined easily. 

5 Renormalization 

This section is devoted to renormalizing the 3 and 4-point contributions to F and G. First, all of the 
3-point contributions must be put in the same form so as to be easily combined. We can then localize 
the ultraviolet divergent pieces, and combine them with the 4-point divergences. By reading off the 
correct counterterm coefiicents we remove all divergences, and are left with the fully renormalized 
structure functions of the vacuum polarization. 

5.1 Converting to Functions of Ax 

The expressions we derived in section 4 contain many indefinite integrals of products of derivatives 
of the three propagator functions iA^^/), Biy) and Giy). Each of these products consists of a few 
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powers of y{x; x') which are singular at coincidence {x'^ = x^) and whose coefficients are nonzero for 
D = 4, plus an infinite series of less and less singular powers of y whose coefficients vanish for D = A. 
Because the vacuum polarization is used inside the 4- dimensional integral of the quantum-corrected 
Maxwell equation (3), the only terms which require dimensional regularization are those which are 
at least as singular as l/y^. Any less singular term can be evaluated for D = 4, at which point most 
of the tedious infinite series contributions vanish. Recalling as well that y{x\ x') = H'^aa'Ax'^, we can 
make the following simplifications: 

For Fbb and Gbb we use the identities 



I[BB] -^ 

I[B^ 

I^[B^ 

I'[B^ 

P[BB] -^ 



r^ (f - 1) 



-)■ 



-)■ 



-)■ 



167r^(aa')^-2Aa;2^-4 

r^ (f - 1) 

327r^(aa')^-'Ax2^-4 

p_2)2 r2(f-i) 

-1) 



{D - 1) 6A7iD[aa 
{D - 2) r2 ( 



+ 



2567r^ 



D 

~2 



laa' 



r\D-2 



Ax^^- 



{D - 1) 647r^(aa')^-2Ax2^-4 
H^ 1 



967r4 aa'Ax2 

D r2( 



1) 



{D - 1) 64TTD{aa')^-^Ax^D-^ 



(106) 
(107) 
(108) 
(109) 
(110) 
(111) 



where we note that for this set of identities only (110) could be put in the D = 4 limit. For Fbc and 
Gbc most of the terms will go to zero and there is only one identity needed 



I[BiB-C))]^ 



(D-4) r2(f-i) 



12^ 



7r^(aa')^-2Aa;2D-4 



(112) 



For Fab and Gab we find that we can take D = 4 in most of the terms 



B{iAA-B) -^ 

I[{iAa - B)B] 

I[{iAa - B)B] 

I\iAa-B)B] 

P[{iAa-B)B] 

I\iAa-B)B] 

kB -^ 



H^ 



327r4 



—?■ 



-)> 



-)■ 



-)■ 



-)> 



i + In l-H'Ax^ 
2 \4 

- + In f-H^Ax^ 
2 V4 



1 



327r4 
(D-2)r2(f-l) 



aa'Ax'^ 

1 
aa'Ax'^ 
1 



Q4tt^ 



[aa 



i\D-2 



Ax 



2D-4 



H^ 



627r4 aa'Ax"^ 



647r4 

_JP_ 

~32^ 



In 



4 
2 + In 
1 



H^Ax^ \ + u 
1 



:H^Ax'' 



1 



aa'Ax"^ 



327r4 aa'Ax"^ 



(113) 

(114) 

(115) 

(116) 
(117) 

(118) 

(119) 
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Notice that (115) is the only identity that needs to be kept in D dimensions. Lastly for Fcb and Gcb 
there is again only one relevant identity since most of the terms will be zero 

(D-4) r2(f-i) 



I[[C-B)B]^ 



128 



n^iaa' 



■\D-2^^2D~A ■ 



(120) 



5.2 Finding the Finite and Divergent Parts of F and G 

It is easiest to renormalize the structure functions term by term; thus we will demonstrate the pro- 
cedure for one term and the reader can extrapolate from the example to derive the rest of the terms. 
First, it is useful to note the necessary identities for renormalization. The terms proportional to 
l/Ax^ are already integrable and do not need to be renormalized. There are also terms proportional 
l/Ax^"^"^ and 1/Ax'^^~^; these will need to be renormalized. We use dimensional regulation to 
partially integrate these terms until they are integrable 

1 d^ 1 



Aa;2-D-2 
1 



A{D - 2)2(D - 3){D - 4) Ax^D- 
d^ 1 



Ax 



2D-4 



2(D-3)(D-4) Aa;2^-6 ' 

but it is clear that these identities contain a divergence in the factors of jj^:^- 
divergence by adding zero in the form 

?47r^/2 



d' 



1 



X 



Ax^ r(f-i) 

Then by adding (123) to the divergent parts of (121) and (122) we find 

d^ 1 i47r^/2 /i^-4,5^(a;-a;') d^ ( I 



(121) 

(122) 
We can localize the 

(123) 



[D - 4) Aa;2S- 



r(f-i) 



^~^6^{x-x') 



/i 



D-4 



(D-4) (Ax2^-6 
Aa;2 



Ax^-2 
0{D-4) 



r(f-l) (Z^-4) 2 

where the factor of fi is added for dimensional consistency. We may now begin our example: 
We consider the eighth term in Fbb 



(124) 



I i/y • i// 



{D^-D- 4) 



{D-2) 



K\v^ + d,d',\ 



Applying identity (HI) we find 

{D^-D 



-^^exl^'*'! X J 



{D-2) 



%\W^ + dod',) 



,aa 



D 



'\D-3t2 



r[BB] 



r(f-i; 



6ATT^{D-l)aa'Ax^^-^ 



(125) 



(126) 



We can now use identities (122) and (124) to break this into its divergent and finite pieces. We find 



-'^ex,divl^'''i X ) 



D{D'^-D-A) K^ fD 

8{D-l){D-2){D-3) IGtt^ V^' 

,4vr^/2 



-,d' + 2-d, 
a^ a 



-'^ex,finitel^''') •'' J 



487r' 



(V^ + d,d',) 



-.9' 



aa 



(D-A) 

ln(;u2Aa;2 

Aa;2 



^^-^r (a; - x') 



(127) 
(128) 
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where in deriving the divergent part we have used the fact that, in conjunction with the delta function, 
dod'o^ — )■ 2— (9o — T<9o- Renormahzing all other terms will follow a very similar procedure. 

Once all of the terms have been renormalized and combined we finally arrive at the full result for 
the structure functions coming from the 3-point diagram, for convenience we split the results into 
their finite and divergent pieces 



i^3,div(a;;x') 



{D^ - QD^ + QD^ + 3^2 + 32D 






{D - 1){D - 2){D - 3){D - A) 



56)^.^ 



D 



8{D-A){D 



a 



-,d' + 2-do 



a 



(129) 




9^ 



'\n{i2'^Ax 


')] 






Ax^ 






H^ 


A 


2 \n{lH^Ax^y 

x"^ Ax'^ 


6 



+ 



3Ax2 



327r2 



1927i^aa' 

\n{j^'^Ax'^ 

Ax2 

In (iH^Ax^) 

A^^ 



3 IndH^Ax"^) Ar]^ 1 1^2, /..2a 2^ 
+ 2 ^\ , ^ + -^do-^ + -V \n{H^Ax^' 

2 Ax^ 4 



\n{jj,'^ Ax'^ 
Aa;2 



Gs^divix; X ) 



GsfinitcyX] X ) 



Ax^ Ax"^ 

[D^ - bD^ + AD^ + 11^2 _ 15/} _ 4) 
8(D - 1){D - 2){D - 3){D - 4) 

ln(/i2Ax2)' 



(130) 




:i3i: 



, 1 

Aa;2 


-¥' 


l^' 


Ax2 [ 



Aa;2 



+ d' 



In (i//2^3,2) 



Ax2 



In(ii72^x2) ^ ^ 



Aa;2 



iv2ln(ij2Ax2^ 



(132) 



In the fiat space limit iJ — )• 0, a = a' = 1 these structure functions recover the old result [12]. We are 
now ready to combine the 3-point and 4-point contributions to find the full vacuum polarization. 

5.3 Our Full Result 

We will now find the appropriate counter term coefficients such as to cancel all divergences, and be left 
with the full renormalized vacuum polarization. Recall that we can actually absorb all of the 4-point 
contribution, minus the terms proportional to ln(a), into the counterterms. Thus to find the full 
counterterm coefficients we simply have to add the divergent coefficients of the 3-point contribution 



28 



to the the 4-point contribution. From (68), (69), and (129) we see 
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from (129) we find 
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and from (70), (71), and (131) we find 
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(135) 



Using these coefficients in the counterterms (52) and (53) we have sufficiently removed all divergences 
and arrive at the final renormalized form for the structure functions 
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These structure functions comprise our main result. 
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6 Hartree Approximation 

Solving the quantum-corrected Maxwell equation (3) is an exercise of comparable difficulty to the one 
we have just concluded, so it will appear in a separate publication. In the meantime, we can gain a 
qualitative understanding of what the result might show by applying the Hartree approximation [48] . 
This has been used previously to study the effect of charged inflationary scalars on photons [49] and 
(especially relevant to the current problem) the effect of inflationary gravitons on massless fermions 
[7]. In each of the previous cases the Hartree approximation gave the correct spacetime dependence 
of the one loop correction to the mode functions and the correct sign relative to the tree order result. 
The Hartree approximation to (3) consists of replacing the Heisenberg operator field equation by 
its expectation value in free graviton vacuum and then expanding in terms of coincident graviton 
propagators. 
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Here JJ^^^P"^P'y^ denotes the tensor factor of the 2-graviton-2-photon vertex, given in expression (19). 
Most of the coincident graviton propagator is a divergent constant; the secular effects for which we 
are searching derive from only the logarithm part of the A- type propagator (25). At this point we 
can also take D = 4,^ so our Hartree approximation to the effective field equation (3) is. 
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Recall that the A- type tensor factor was defined in expression (32). 
Substituting (19) and (32) in (141) gives a simple result. 
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(142) 



Here a barred index on any tensor indicates that its 0-component has been suppressed, for example, 
V^ = ri^^Vy = V^ — 5qV^ . We can distinguish in expression (142) between the cases of z/ = and 
u = i. The constraint equation is effectively multiplied by a secular factor. 
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^One might worry about factors of ln(a) which could arise when a divergent constant multiphes a^^^. However, 
we can see from expressions (52) and (53) that the very same factor of a^~^ multiphes the counterterms which absorb 
divergences from the 4-point diagram. So there can be no secular contributions from this source and we may as well 
drop the divergent constants and take D — 4. 



30 



This has no effect on dynamical photons ahhough it would lead to a secular screening of a point 
charge of the sort reported by Kitamoto and Kitazawa [51] provided there is no compensating secular 
factor on the charge density. The equation of relevance for dynamical photons is u = i, 
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The one loop correction to the effective field equation of course fixes only the one loop corrections 
to the field strength. We therefore expand in powers of k^, 
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Equations (144) and (145) imply the following relations for the one loop field strengths, 
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With the f/(l) Bianchi identity, the leading secular behavior is. 
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We see that the one loop correction to the electric field strength of a photon tends to cancel its 
classical value whereas the one loop correction to the magnetic field strength dies off. 

7 Discussion 

We have used dimensional regularization to compute the one loop quantum gravitational contribution 
to the vacuum polarization on de Sitter background. We first calculated the 4-point contribution in 
section 3, and then derived the much more cumbersome 3-point contribution in section 4. Each result 
was expressed in the form (40) as the sum of two transverse projection operators acting on structure 
functions. In sub-section 2.5 the relevant BPHZ counterterms (47) were also reduced to this form, 
resulting in expressions (52-53). Renormalization was implemented in section 5 to give our final 
results (136) and (137) for the structure functions F{x;x') and G{x]x'). 

Our ultimate goal is to study how infiationary gravitons affect electrodynamics using the quantum- 
corrected Maxwell equation (3). Specializing to de Sitter in conformal coordinates, substituting our 
form (40) for representing the vacuum polarization, and partially integrating the primed derivatives, 
allows us to express the quantum-corrected Maxwell equations as. 



d^F'^^i^x) + d^ (Tx'l iF{x;x')F''^'{x')+iG{x;x')F^^{x') 



J^ . 



(151) 



(Recall that a barred index is purely spatial.) Equation (151) can be employed the same way one 
uses the classical Maxwell equation to study dynamical photons (J^ = solutions) and the electric 
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and magnetic fields induced by standard sources. We have already done this for the one vacuum 
polarization from gravitons on flat space background [12]. Closely related studies have also been 
made of the effects that inflationary scalars have on dynamical photons [3] and on electrodynamic 
forces [4]. 

The actual implementation of this program requires solving integro-differential equations in the 
context of the Schwinger-Keldysh formalism [50]. That is an project comparable to the one we have 
just completed, so it will be deferred to a separate work. However, a simple estimate of what it might 
give was derived in section 6 by making the Hartree approximation [48] to localize the effective field 
equation. We find that the one loop electric field strength (149) of dynamical photons experiences 
a secular growth which tends to cancel its free field value, whereas the one loop correction to the 
magnetic field strength (150) dies away compared to its free field value. It is interesting to note that 
the magnetic response to inflationary scalars also seems to be subdominant to the electric repsonse 

[4]. 

Working out what the full equation (151) gives for dynamical photons is important to check the 
observation in [9] that the spin-spin interaction between gravitons and fermions seems to explain 
why inflationary gravitons cause the fermion mode function to grow [7] whereas they have no secular 
effect on the mode function of a massless, minimally coupled scalar [8]. Another important exercise 
is to work out the effect of inflationary gravitons on the electric field of a point charge. This is the 
natural way to check the surprising claim of Kitamoto and Kitazawa that infrared gravitons screen 
sub-horizon interactions during inflation [51]. 

Before closing, we should comment on the gauge issue. The vacuum polarization requires fixing 
both the ^(1) and diffeomorphism symmetries, and the manner in which this is accomplished can 
affect the result. Our previous study of gravitons on flat background revealed no dependence upon 
the choice of electromagnetic gauge, but a huge variation with the gravitational gauge [12]. We 
believe there is not likely to be any gauge dependence in the leading secular infrared effects one finds 
from de Sitter gravitons because the spin two part of the graviton propagator has the same infrared 
logarithm term in any gauge [30, 28]. Note that it is perfectly possible for a IPI function such as the 
vacuum polarization to change with the gauge, while a particular feature of its dependence on space 
and time is the same in all gauges [52]. That is precisely what happens with the pole terms of IPI 
functions in flat space quantum field theory, and we suspect that the same applies for the leading 
secular dependence on de Sitter. 
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